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Abstract. We give an overview of the two different methods that have been introduced in 
order to describe the dynamics of constrained quantum systems; the symplectic formulation 
and the metric formulation. The symplectic method extends the work of Dirac on constrained 
classical systems to quantum systems, whereas the metric approach is purely a quantum 
mechanical method having no immediate classical counterpart. Two examples are provided 
that illustrate the nonlinear motion induced by the constraint. 



1. Introduction 

A question that has recently been asked is how does one implement unitary motion on a quantum 
system that is subject to a set of constraints? In classical mechanics, constrained evolutions 
are often constructed using an approach developed by Dirac [TJ [2], which uses the symplectic 
geometry of the classical phase space. Recently, several authors [31 HI [5] have worked on extending 
the classical theory of constraint to quantum systems, using the phase space formulation of 
quantum mechanics. There are, however, pros and cons in this approach: It is for example only 
applicable to systems subject to an even number of constraints. To address this and other issues, 
an alternative approach for treating constrained quantum systems was introduced in [6] . In this 
paper we review both these approaches and explore through examples the kind of nonlinear 
motion that results from constraining a unitary evolution in quantum mechanics. 

The classical method uses the property that the classical phase space is a symplectic manifold 
endowed with a symplectic structure. The quantum mechanical phase space has the geometry 
of a Kahler manifold, with both a symplectic structure and a metric structure. As the method 
considered by [SHUE] uses the symplectic geometry of the quantum phase space, it was given the 
name the symplectic approach to quantum constraints [5] • The alternative formalism introduced 
in [6] makes a novel use of the metric structure of the quantum state space, thus referred to as 
the metric approach to quantum constraints. 

In this paper we will first review some of the main features of the phase space formulation of 
quantum mechanics. We will then highlight the key features of both the symplectic and metric 
approaches to implementing constrained unitary motion and discuss the resulting evolutions by 
looking at two examples. 



2. The phase-space formulation of quantum mechanics 

The phase space description of quantum mechanics is based on a projective Hilbert space 
formulation (see for example El H2 QUI [HI Q2J H3 O [15] and references cited therein) . The 
idea is as follows. The expectation value of any observable in quantum mechanics is independent 



of the norm of the state vector 



x) ~ X\x) 



A eC-{0}. 
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A state in quantum mechanics can hence be viewed as a ray through the origin of the Hilbert 
space TL n+1 , and we can, without loss of any physical information about the system, consider 
quantum mechanics in the space of rays. The space of rays correspond to the complex projective 
space P ra , where each state of the system is now given by a point. In real terms P n is an even 
dimensional manifold T 2n equipped with a Riemannian structure given by the Fubini-Study 
metric g a b and a compatible symplectic structure f2 a 6. The transition probability between two 
states in T is given by the associated geodesic distance between the two points. It is interesting to 
note that the probabilistic aspects of quantum theory are hence associated with the underlying 
metric geometry of the state space. In comparison the dynamical aspects of the theory are 
captured by underlying symplectic geometry. Physical observables are functions on T given 
by the expectation value of the corresponding Hilbert space operator at each point x. The 
Hamiltonian, for example, is given by 



and the corresponding eigenvalues are the values of H(x) at each fixed point. Unitary evolution 
is governed by the Schrodinger equation, which takes the form of Hamilton's equations: 



where O is the inverse symplectic structure given by Vt^O* = —5^. Here, and throughout the 
paper, there is an implied summation over repeated indices. We parameterise the state space 
by the "action-angle" parameterisation (cf. [16j) so that the equations of motion (jl]) take a 
particularly simple and familiar form. If we let the state \x) in energy basis be given by 





The eigenstates of the Hamiltonian are given by the fixed points of H: 
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^ = n ab v b H( x ) 
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then the expectation of the Hamiltonian operator H = ^ - Ej\Ej){Ej\ is given by 





(6) 



The equations of motion then take the form 




dH(qi,pj) 
dqi 



(7) 



with solutions qi{t) = qi(0) + u>it and Pi(t) = Pi(0). 
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Figure 1. The constraints = define a 
subspace of T. The evolution of the system can 
be constrained to this subspace using either the 
symplectic or the metric method. When using 
the symplectic approach the modified equations 
of motion take the same form as the original one, 
but with a modified symplectic structure £l a b- In 
the metric approach the system is constrained by 
removing all components orthogonal to the surface 
from the vector field x a = Q, ab VbH. 



3. Constrained dynamics 

Let us now assume that the evolution of the quantum system is subject to one or more 
constraints. The type of constraints that we will be considering can be written in the form 



<S> l (x) 



0, 



1,2,... ,N, 



(8) 



where N is the total number of constraints. Equation ([8]) defines a subspace <3? of T, given by 
the intersection of each of the N subspaces defined by = 0, onto which the evolution will be 
constrained, see figure [TJ 

We shall here consider two different types of constraints. The first type is given by the 
conservation of a set of observables, i.e. 



{x\$ l \x) 
(x\x) 



(9) 



where is the Hilbert space operator corresponding to the observable $\ The second type 
of constraint is when the evolution of the system is constrained to an algebraically well defined 
subspace of V that does not correspond to the conservation of any observable. An example of 
this is to constrain the motion of a pair of spin-^ particles that are initially disentangled to 
remain disentangled. There is no observable that corresponds to this property of the system, 
but as disentangled states form a subspace Q = P 1 x P 1 [151 CL7] of the overall state space P 3 
(see figure [2]), we can constrain the system algebraically to remain on this subspace. We shall 




Figure 2. The state space of 2 spin-i particles. 
The disentangled states form a subspace of the 
total state space given by Q = P 1 x P 1 , where P 1 is 
the complex projective line. Further details of the 
state space structure of this system can be found 
in [13 IT]. 



return to this example and look at it in more detail in the next section where we give an outline 
of the symplectic approach to implementing quantum constraints. 



3.1. Symplectic approach 

Let us first consider the symplectic approach and follow the methodology introduced in [5]. In 
order to implement the constraints ([8]) the equations of motion are modified such that 

— = n ab v b H(x) + \in ab v b &(x), (io) 

where Aj are Lagrange multipliers associated with the constraints <& l (x). In some cases these 
Lagrange multipliers can be found explicitly by considering <I> 1 = 0. By expanding <3? 1 using the 
chain rule <I> 1 = x a V a <3? J we find that the Lagrange multiplier is given by: 

\i = u ji n ab v a &VbH, (ii) 

where wy, assuming it exists, is the inverse of 

Jj ._ Qab Va $i Vb $j ( 12 ) 

Note that when the constraints are given by conserved observables ©, then uf 3 correspond to 
the commutator between the two observables. By substituting the Lagrange multipliers back 
into (I10p we find that the equations of motion can be rewritten in the form [5] 

^ = n ab V b H(x), (13) 

where Q ab is the induced symplectic structure on the constraint subspace & l = given by 

n ab = n ab + n^n^Vc^Vd**'. (14) 

In other words, the equations of motion for the constrained system take the same form as 
for the unconstrained system. Only the symplectic structure is no longer the globally defined 
f2 a b, but the symplectic structure induced locally on the subspace <!>* = (see figure HJ. A 
consequence of this is that we can no longer expect the resulting motion to be unitary, as the 
induced symplectic structure VL ab has a nonlinear term. In order to illustrate the kind of motion 
that could result from this kind of constraint we consider the example of two spin-^ particles 
constrained to remain disentangled, as mentioned above. Let us first note however that the 
symplectic approach is not applicable to all systems. In order to use it, one requirement is that 
N, the total number of constraints being imposed, must be even. This follows by construction as 
otherwise the symplectic geometry is not preserved on the constraint surface. Another condition 
that must be satisfied in order to find the lagrange multipliers is that the inverse of (I12p must 
exist. 

Example 1: Let us now consider the example of two spin-^ particles with Hamiltonian: 

H = -Jd 1 ®& 2 - B(&l®l 2 + li®a%), (15) 

where the subscripts 1 and 2 labels the two particles, J is the strength of the spin-interaction, 
B is the strength of the external magnetic field (orientated along the z-axis) , and a z is the Pauli 
spin matrix in the z-direction. We further require the initial state to be disentangled, and we 
wish to constrain the system so that it remains so. In other words we want the state of the 
system to remain on the surface Q = P 1 x P 1 as it evolves. To find the equations of motion we 
follow the method outlined above, first finding (fT2j) . taking its inverse, using it to evaluate the 
Lagrange multipliers and hence finding Q ab (please see [5] for full details of this calculation). 

Since the complex projective line is isomorphic to the two-sphere P 1 ~ S 2 , we can view Q as 
the product space of a pair of 2-spheres Sf x Sf,- Changing coordinates into spherical coordinates 




Figure 3. Four examples of 'snapshots' of the vector field and its integral curves as two spin-^ 
particles are constrained to remain disentangled so that their evolution lies on the product space 
S 2 x Sjj. The figures show the vector fields on S 2 as the second particle evolves through the 
states (0 2 = § ,0 2 = f) (top left), (0 2 = f ,02 = 0) (top right), (0 2 = if, 02 = if) (bottom left) 
and (0 2 = ^-,fo = §) (bottom right), on S|. 



(see [5]) we can visualise the evolution on the surface of the 2-spheres. The resulting equations 
of motion consist of four coupled non- linear differential equations: 

01 = sin(0i - 2 ) sin0 2 [(wi - u 2 ) cos0i + u 2 — u 3 ] 

2 = sin(0i - 4> 2 ) sin 0i [(w 2 - wi) cos 6 2 - uj 2 + u 3 ] 

01 = |[— + (cj 2 - ^) cos 2 + (|wi - uj 2 - 2^3) cos 0i + [cos (01 

x (2(cj 3 — uj 2 ) sin 2 0i cos 02 + (uj\ — uj 2 )(cos 2 B\ — cos 2 02))] 

02 = + (uj 2 - ^) cos 0i + (§wi - uj 2 - 2uj 3 ) cos 9 2 + [cos(0i 

x (2(cj 3 — uj 2 ) cos 0i sin 2 2 — {ui\ — uj 2 )(cos 2 B\ — cos 2 02))] 

It is difficult to solve these equations of motion in closed form. We will therefore explore 
the solution numerically by looking at various 'snapshots' of the vector field and associated 
trajectories as the system passes through certain points. In figure we see four examples of 
this. The top left picture shows us the vector field on S 2 as the evolution on S 2 passes through 
the point 2 = 02 = \k. As the overall state of the system changes, we evolve on both spheres 
and as we change state on S 2 the vector field on §| changes and vice versa. To fully explore 
the solution numerically, we require an 'interactive simulation'; figure [3] merely gives a flavour 
for the often complicated trajectories that result from constraining the system to remain on 
Q. The strength of the formalism is that the solutions are nonetheless fully tractable, at least 
numerically. 



- 02 ) /(sin 0i sin 2 )] ^ 

- 02 ) /(sin 0i sin0 2 )] 



3.2. Metric approach 

In the previous section we found that there is a way to implement quantum constraints using 
the underlying symplectic geometry of the phase space. Let us now look at how some of the 
limitations of the symplectic formalism can be overcome by using the metric geometry of the 
phase space. 

Another way of constraining the evolution of a system whose initial state lies on the surface 
<& = is to remove all the components of the vector field normal to it [B] (see figure [T]). This 
will force the system to remain on $ = 0. To implement this constraint we still use the method 
of Lagrange multipliers, but we substitute the symplectic structure in the second term of (fTU|) 
with the metric structure. Hence the modified equations of motion become: 

c\r a 

= n ah V b H-\ ig ab V b <5>\ (17) 

at 

where 

g ab ig 

the inverse of the Fubini-Study metric. Because we are now making use of the metric 
geometry of the constraint subspace, we no longer require that the total number of constraints 
N be even. The Lagrange multipliers Aj can again be found explicitly in a similar way to the 
previous approach by considering <E> 1 = 0. Following the method in [6], we hence find that 

Xi = M ij SV*V a &V b H, (18) 

where is the inverse of the matrix 

M ij ^v o $ i V 6 * J '. (19) 

For the case when the constraints correspond to a family of conserved quantum observables, 
the matrix M tJ corresponds to the anticommutators between the observables. We hence require 
det(M lJ ) ^ in order to find the equations of motion: 

x a = Q ab V b H - M ij n cd V c &V d Hg ab V b <I> i . (20) 

It is no longer in general the case that these can be rewritten in the form x a = Cl ab V b H, where 
Q ab is the inverse of a modified symplectic structure. We would, however, expect that systems 
that can be treated with both methods return the same result. It turns out that this is very hard 
to verify in general. In [6], a sufficient condition was derived under which the metric approach 
reduces to the symplectic approach with the modified equations of motion taking the form of 
(|13p . (see [6] for further details on this). We therefore take the view that the metric approach is 
the more general way of treating quantum constraints, but under certain conditions it reduces 
to an extension of the classical method to quantum systems as in the symplectic approach. We 
shall now look at one of the simplest examples of a system subject to only one constraint. 

Example 2: Let us consider the example of a single spin-^ particle with Hamiltonian H = a z , 
where a z is a Pauli spin matrix, constrained such that 

*(*) = ^ (21) 

is conserved [6] . This example is interesting as we have an observable that does not commute with 
the Hamiltonian. Following the method outlined above, the equations of motion are obtained 
by considering (|19|) . taking its inverse, finding the Lagrange multipliers (|18|) and hence the 
equations of motion (|20p . Just as in the previous example we look at the problem in spherical 
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Figure 4. The vector field resulting from a single spin-^ particle with Hamiltonian H = a z 
being constrained such that A[x) = ^^j^ is conserved [6]. The yellow lines are integral curves 
of the vector field and the red lines correspond to fixed points. 



coordinates to more easily visualise the resulting evolution as trajectories on the surface of the 
2-sphere. The modified equations of motion are given by [6] 

sin(20) sin(2<£) \ , • 2 cos 2 9 cos 2 <j) , . 

t, — , and, <p = ^ — . (22) 

1 — sin 9 cos 2 (j> J 1 — sin 9 cos 2 <f> 

Figure S] shows the resulting vector field and some of its integral curves. We recall that unitary 
unconstrained motion for this system corresponds to rigid rotation around the z-axis. As we 
are now forcing the system to remain on (a x ) = 0, corresponding to circles around the x-axis, 
the resulting solution has a very interesting fixed point structure. Instead of only having two 
fixed points, the poles along the z-axis, as in the unconstrained case, the state space is now 
effectively partitioned into quarters of the sphere that do not interact with each other. Systems 
whose initial state is not a fixed point hence evolve along the half circular trajectories towards 
the fixed points along the equator. 
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